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Abstract: 

We present calciilations of axial-current matrix elements between various heavy-meson and heavy-baryon states to 
the next-to-leading order in heavy hadron chiral perturbation theory in the p-regimc. When compared with data 
from lattice computations or experiments, these results can be used to determine the axial couplings in the chiral 
Lagrangian. Our calculation is performed in partially-quenched chiral perturbation theory for both SU(4|2) and 
SU(6|3). We incorporate finite-size effects arising from a single Goldstone meson wrapping around the spatial volume. 
Results for full QCD with two and three flavours can be obtained straightforwardly by taking the sea-quark masses to 
be equal to the valence-quark masses. To illustrate the impact of our chiral perturbation theory calculation on lattice 
computations, we analyse the SU(2) full QCD results in detail. We also study one-loop contributions relevant to the 
heavy hadron strong-decay amplitudes involving final- state Goldstone bosons, and demonstrate that the quark-mass 
dependence of these amplitudes can be significantly different from that of the axial current matrix elements containing 
only single hadron external states. 
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I. INTRODUCTION 



The physics of b hadrons is an important and active field of research, both experimentally and theoretically. B 
mesons have played an important role in our understanding of flavour physics in the Standard Model (SM) and 
its possible extension. The on-going LHCb experiment and possible future B factories will produce significantly 
improved experimental information for B mesons which will, in turn, lead to better constraints on the relevant SM 
parameters or reveal deviations from the SM. In addition, a large amount of polarised single-bottom baryon data will 
be produced. This will allow extensive studies of the spectrum and the decays of these baryons. Since the baryons 
carry different spin quantum numbers, they may offer additional opportunities for probing the coupling structure of 
physics beyond the SM. In performing such investigations, it is necessary to compare experimental results to precise 
theoretical calculations in which non-perturbative strong-interaction effects are well controlled. This is becoming 
achievable because of the progress in Lattice QCD. 

Calculations in Lattice QCD are often performed at unphysical light-quark masses due to the limited computing 
resources. In order to obtain high-precision theoretical predictions for spectral quantities and matrix elements, it is 
essential to use chiral perturbation theory (xPT) to extrapolate to the physical quark masses. For systems of hadrons 
containing a single valence 6 or 6 quark, the relevant chiral effective field theory is heavy-hadron xPT (HHxPT) [I]- 
[5]. In addition to the low-energy constants in the chiral Lagrangian of the Goldstone boson sector, there are three 
unknown coupling constants in this effective theory at the leading order (LO). These constants, defined explicitly as 
5i,2,3 in Eq. (24 1 in Secjllj accompany axial couphngs of heavy hadrons to the Goldstone boson sector and appear in 



all chiral extrapolations using HHxPT. Therefore, the accurate determination of 31,2.3 is one of the most important 
tasks in the Lattice QCD calculations for 6-physics phenomenology. 

In this work, we compute the matrix elements of the quark-level axial currents, 

Jud.fi = dlfi^bU, and Jus,^ = sjfj^j^u, (1) 

between various heavy-light meson and single— 6 baryon states to the next-to-leading order (NLO) in HHxPT. In 
particular, we calculate the relevant one-loop contributions to these matrix elements. When compared with data from 
lattice calculations or experiments, our results can be used to extract the above-mentioned three axial couplings in 
HHxPT. Our calculation is performed in partially-quenched chiral perturbation theory (FQ^PT) using the super- 
symmetric formulation 6 , for both SU(4|2) and SU(6|3). The "fuU-QCD" limit can be taken straightforwardly from 
our results by setting the sea-quark masses to be equal to the valence-quark masses. Our one-loop computation is 
carried out for finite spatial volume in the p- regime ^, following the same method as in Refs. [HIID]- As pointed out in 
Ref. (llj . in heavy-light meson systems, finite- volume effects arising from higher-order terms in the chiral expansion 
can be estimated. This requires high-precision information on the B*—B—t: coupling beyond that which is currently 
available. Nevertheless, such higher-order effects are insignificant for current and future lattice calculations, since 
computations with small pion masses in large volumes are becoming standard. 

In this paper, we present our results in the isospin limit. However, in the case of SU(6|3), we include the SU(3) 
breaking^ effects, both in the external states and in the axial currents. At NLO in HHxFT, the axial-current matrix 
elements for heavy hadrons can be written in the general form 

g{l + g^L + cpL' + L") + analytic terms, (2) 



where g and g' are variously gi, g2 and 53 in Eq. ([24]), and L, V and L" are the contributions from one- loop diagrams. 
The determination of g\ using lattice QCD has been attempted by various groups [T^HIZ] . However, the correct quark 
mass dependence (based on the symmetries of QCD) of the axial matrix elements was not previously known. Using 
the current work, extrapolations to the physical quark masses can be made rigorously. 



This paper is organised in the following way. Section [TI] contains an introduction to HH^PT. In Sec. |III[ we first 
present the general structure of the one-loop contributions to the axial-current matrix elements, before giving the 



^ Studies of the heavy-meson systems in the e-regime can be found in Ref. [71[8]. 

^ More precisely, we consider identical SU(3) breaking effects in the sea, valence and ghost sectors of SU(6|3) but for simplicity, refer to 
this as SU(3) breaking. 
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results in the case of SU(2) in Sec.|lV| Results for SU(4|2) and SU(6|3) HH^PT are presented in Sec.|v) emphasising 



the quark flavour flow picture. In Sec. 
also computed before we conclude. Technical details of the results are included in the appendices. 



VI the strong-decay amplitudes involving flnal-state Goldstone mesons are 



II. HEAVY HADRON CHIRAL PERTURBATION THEORY 



The partially quenched (PQ) chiral Lagrangian^ for the Goldstone mesons is 

£g - ^str [(a^Et) (d^^) + + x^S] + [a{d^<S>,){d^<fo) - M^cf^j , (3) 

where S — exp(2i$//) is the non-linear Goldstone particle field, with $ being the matrix containing the standard 
Goldstone fields in the quark-flavour basis. We use / — 132 MeV. In this work, we follow the supersymmetric 
formulation of PQ chiral perturbation theory (PQxPT) [6 . Therefore under SU(4|2)l d) SU(4|2)r or SU(6|3)l d 
SU(6|3)r, E transforms as 

S (Jl S Ul (4) 

where 



Ul e SU(4|2)l or SU(6|3)l, 



C/r e SU(4|2)r or SU(6|3)r. (5) 
The symbol "str" in the above equation means "supertrace" . The variable x is defined as 

X = 2BoMg, (6) 



where Bq is a low energy constant related to the chiral condensate and, in the isospin limit, the quark mass matrix, 
Mq is 

Mq = diag ( m„^m,„ m„/ ^to„^ rnu^ rriu) , (7) 

valence sea ghost 

in the SU(4|2) theory, and is 

Mq = diag( mM , to„ , , m^/ ,mu', rris'^ , m„ , , rus ) , (8) 

valence sea ghost 

in the SU(6|3) theory. We keep the strange quark mass different from that of the up and down quarks in the valence, 
sea and ghost sectors. Notice that the flavour singlet state $0 = str($)/-\/6 is rendered heavy by the U{1)a anomaly 
in PQQCD [191 [20] and can be integrated out, resulting in residual "hairpin" structures. 



The inclusion of the heavy- light mesons in chiral perturbation theory was first proposed in Refs. [HEj: with the 
generalisation to quenched and partially quenched theories given in Ref . [HI [52] . The 1 /Mp and chiral corrections 
were studied by Boyd and Grinstein [23j. The B and B* meson fields appear in this effective theory through the 
"superfield" 



= K^7^-S,75) V^, (9) 



^ In this paper, we only address situations where there are no multi-particle thresholds involved in loops. Therefore, in spite of the sickness 
pointed out in Ref. jl8| . we can still use the Minkowski formalism of PQ chiral perturbation theory. 
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where t;^ is the 4- velocity of the meson fields, Bi and B* ^ annihilate pseudoscalar and vector mesons containing 
an anti-6 quark^ and a light quark of flavour i. Under the heavy-quark spin transformation S^, and the unbroken 
light-flavour transformation U{x), the field H^''^ transforms as 



iff (x)^C//"(x) Hf\x)S^- 



(10) 



Also, the conjugate field, which creates heavy-light mesons containing an anti-6 quark and a light quark of flavour i, 
is defined as 



(b) 



and transforms under Sh and U (x) as 



(x) {U^)\{x). 



(11) 



(12) 



The introduction of the single— 5 baryons to xPT was pioneered by authors of Refs. [3HS], and the effective theory 
was generalised to the PQ scenario in Ref. |25j . Since the two valence light quarks in such baryons may carry total 
spin quantum numbers^ s; = or = 1, there are two types of heavy baryons. At the quark level, these two types 
of baryons carrying light flavours i and j are described by the interpolating flelds 



a, a B.b , Q.b a, a 

q, q'-' +q'-' q. 



iabc {Cl5)ap foi' «i = 0: 



57. 



Q.b Q,a 
1, 9^ 



Eafjc {.CY)aB = 1' 



(13) 



where C is the charge-conjugation matrix, a, /3 and 7 are the Dirac indices and a, b and c are colour indices. In full 
QCD, the T flelds are anti-symmetric and the S flelds are symmetric under the exchange of the light flavour indices. 
In the PQ theory, the flavour structure of these interpolating flelds has the properties 



Tij — ( ^Y'^^^Tji, 



where 



1 when i € valence and sea, 
when i S ghost. 



(14) 



(15) 



accounts for different statistics of quarks in PQQCD. These flelds transform as 39- and 42-plets under the SU(6|3) 
flavour rotation, while they transform as 17- and 19-plets under the SU(4|2) flavour rotation. The baryon fields are 
included in heavy-hadron chiral perturbation theory (HH^PT) according to the fiavour properties in Eq. (14). In the 
case of Nf = 3 {Nf refers to the number of sea-quark flavours), the pure valence- valence sector of the s/ = baryons 

i 1 /2 

is related to the physical states of A;, and Sj^ ' via 



T, 



(valence— valence ) 




(16) 



where the superscript indicates the 3-component of the isospin. Since the light-light di-quark is of spin-1 in the S'j^j 
fields, such baryons can be in spin 1/2 or 3/2 states which are degenerate in the heavy quark limit. Therefore they 
are best described by the "superfield" 



^ij - 



{v^+r)l5B,,+B* 



(17) 



* We follow the standard notation |24| for the flavour content of B mesons, so that e.g. = _B+ = uh. 

^ The total spin of the light degrees of freedom is a conserved quantum number because of the heavy quark symmetry. 



5 



where Bij and B*^ are spin- 1/2 and 3/2 baryons. In the pure valence- valence sector, 



B. 



(valence— valence 



) = 



1 Y,o \ 



/ y+i _ _ 

VO V^l J_^'-l/2 

V2^b ^-b 



\ ^/2-b ^^b 



(18) 



and similarly for the B*^ fields. The S^J and Tij fields have the same property as ijf'''' under the heavy-quark spin 
transformation Sh- For the unbroken light-fiavour transformation, 



T,j{x) U,\x)U/ix)Tkiix), 



(19) 



with the flavour indices satisfying Eq. (14 1. These S^j and Tij are annihilation field operators and we denote the 
corresponding creation fields by S^j and Tij. 

The Goldstone mesons couple to the above heavy meson and baryon fields in the HHxPT Lagrangian via the non-linear 
realisation 

e = e**/^ = VS, (20) 

which transforms as 

ax) Ul ax) U\x) = U{x) ax) Ul (21) 
The ^ field can be used to construct vector and axial-vector fields 



(22) 



The vector field can then serve as the gauge field in defining the chiral covariant derivative which acts on the heavy 
hadrons, 

P^i?® = df'Hf^ + {¥'')/ Hp, 
OT,, = 9'^7^y + (F'^),'Tfe, + (-1)"'(''^+''*)(1/^)/t,,, 

The leading-order HHxPT Lagrangian is then 



(23) 



C 



(LO) 
HHxPT 



-ztrD 
- 5i tro 



^92 e^..p(5'^«'^A^^'')j + V2g3 [{TA^S^)^ + {S^A^^T)^] , (24) 



where is the velocity of the heavy hadrons, trD[ ] means taking the trace in Dirac space, and ( )j. is the implemen- 
tation of the PQ-theory fiavour contraction rules pS] 

(TFT) J = T^%%j, 



(TY>'S, 



(25) 
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The parameter A^^^^ is the mass difference between the S and T fields with the same light flavour indices, 

A^-^) = Ms- ■ - Mt- .. 

It is grouped together with the definition of other mass parameters in Eq. (Al) in Appendix [A| This mass difference 
is of ©(Aqcd), and does not vanish either in the chiral limit or in the heavy-quark limit. 



The LO Lagrangian for HHxPT contains terms of 0{p) and no light-quark mass dependence. To generate the flavour 
SU(3) breaking effects in heavy-meson and baryon spectrum, which give rise to the mass differences 5^ ^ 
in Eq. (Al), one introduces 

Ax) 



and 8^1, 



C 



HHxPT 



= Al trD 



H, X^Hj 



+ A2 trD 



str {xi) 



+ A3 (5^X4 5^) f H- \i{s^s^)^siT (x^) + A5(rx5r)j + k{tt)^ str (xe) , 



where 



X€ = exe + e^xf^ 



(26) 



(27) 



In the computation of the axial-current matrix elements, the flavour breaking effects in Eq. (261 are formally sub- 
leading compared to those encoded in the pure Goldstone Lagrangian, Eq. (|3]). Nevertheless, we keep them in our 
calculation as they can be numerically significant. 



In this work, we also include the heavy quark spin symmetry breaking term 



X2 
M, 



-trD 



.tin 



(28) 



where Mb is the B meson mass. This counterterm leads to the mass difference between the B* and B mesons with 
the same light flavour. 



M 



which vanishes in the heavy-quark limit. This mass difference is also grouped together with other mass parameters 
in Eq. (Al) in Appendix |A] In principle, there are also such heavy quark spin breaking terms in the baryon sector, 
resulting in mass differences between Bij and B*^ baryons in Eq. ( 17). However, these mass differences are numerically 
much smaller than A^*^) IMI. 



III. AXIAL CURRENT MATRIX ELEMENTS AT THE NEXT-TO-LEADING ORDER 



Applying the Noether theorem to the chiral Lagrangian in the previous section, one can derive the leading-order axial 
currents corresponding to their quark- level counterparts in Eq. ([T]). For matrix elements involving external states of 
single heavy hadrons, the relevant LO currents are 



^^=51 trD 



Cb) f^{+) 



ki 



where the subscript ij means the current changes the light quark flavours from i to j, and 



with the matrices Tij defined as 



(29) 



(30) 



(31) 



where k and I run through all the light-quark flavours in PQQCD. The superscript Nf is the number of sea quark 
flavours and Nf = 2,3 represent the cases of SU(4|2) and SU(6|3) respectively. These leading-order axial currents 



7 



\ 



(a) 



(b) 



(c) 



FIG. 1: One-loop diagrams contributing to the matrix elements of axial currents between heavy hadrons. The dashed lines 
are the Goldstone meson propagators, including the possible "hairpin" structures. The single solid lines denote generically the 
external heavy hadrons, while the double solid lines are the internal heavy hadrons. T hey can be B, B* mesons or Tij, 5, 

(N f) 

baryons. The circled crosses are the insertions of the LO axial current ' given in Eq 



the strong chiral Lagrangian in Eq. ( 24 1 . Diagram (a) is the self-energy of the heavy hadron and it leads to the wavefunction 



29 1, while the other vertices are from 



renormalisation contribution to the matrix elements. Diagrams (b) and (c) are the "tadpole" and "sunset" types respectively. 



generate the LO terms, as well as the NLO contributions via one-loop corrections, in the matrix elements studied in 
this work^. 

There is a significant increase in the number of terms in the next-to-leading order axial currents in HHxPT, and we 
postpone the detailed investigation of these NLO currents to Subsection |IIIB| below. Here we first write down the 
generic form of the chiral expansion of the axial-current matrix elements to the NLO, 



1 



1 



'ij 



{Nf) 
Hi 



(Nf) 



{Nf 



(32) 



where the equality symbol means the matching between (PQ)QCD and the chiral effective theory, /i is the Lorentz 
index, and Jij^^i are the quark- level currents given in Eq. (nj). The fiavour indices are denoted by i and j which are 
not summed in the above expression, and Hi is a heavy hadron state (meson or baryon) containing the light fiavour 



i. The symbols 7^^ 



(Nj) 



(Nf) 
Hi 



and Q^^l^. are results from the tadpole, wavefunction-renormalisation and sunset 



diagrams at one- loop, as depicted in Fig[T] where single and double solid lines represent generically the external and 
internal heavy hadrons while the dashed lines are the Goldstone propagators. The circled crosses in Fig [l] are the 



insertions of the LO axial current Jif^^ given in Eq. (29). The tadpole contributions T}^'^^ are independent of the 
external states, since they emerge completely from the tlavour structure of the currents. 



In Eq. ( 32 1 , we have written the NLO analytic terms as the LO matrix elements times J^^l^j^ ■ In Subsection 
below, we will study these NLO analytic terms, and show that they can be presented in this manner. In this section. 



IIIB 



we examine the analytic terms (polynomials in the Goldstone masses) in the matrix elements in Eq. ( 32 ) for various 
external states. These are encoded in 



and M', 



(Nf) 



Hi^Hi 



The non-analytic contributions arising from the one- loop diagrams will be discussed in Sections |IV| and |V] 



A. Leading-order matrix elements 



Lattice computations are often performed using the baryon interpolating fields in Eq. ( 13 ). Therefore we carry out the 
XPT calculation for the Tij and 5*^^ external states. From our results, it is straightforward to obtain matrix elements 



In addition to the terms in Eq. \29\ , there are other operators in the LO currents arising from the chiral Lagrangian introduced in the 
previous section. Nevertheless, these terms do not appear in the matrix elements to the order we work at. That is, their contributions 
to the one-loop corrections vanish. 
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for physical external baryon states using Eqs. (16) and (18). The leading-order HHxPT predictions for the matrix 
elements studied in this work are 



{^d\Jud!ti \Bu)lo - {B*\jil]f_,\Bu)LO 
{Sdd\jl'^!^\Tdu)LO - V2{S,d\ji%\T,u)LO - V2{Sds\ji%\Tdu)LO = {SsAJL%\Tsu)lo 

{Sdd\jl^'l^\Sdu)-LO = V2{Ssd\Julf,\Ssu)-LO = V2{Sds\j!£lfj_\Sdu)-LO = ('S'ssI >/£^^|5's«)lo 



= -2 5i<, 
= -53 Uf,U, 



'V2 



92 e.^.p U"UP, 
(33) 



where is the polarisation vector of the B* meson, U is the Dirac spinor of the T baryon, and the C/^'s are the 
"superfield spinors" of the 5 baryons. The basis polarisation vectors and spinors satisfy the spin sums 



s = l 

j2u{v,s)U{v^s) = 1±^, 

s=l 
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(34) 



Note that is not a Rarita-Schwinger spinor; instead it contains the degrees of freedom of both the spin- 1/2 and 
spin-3/2 components of the superfield. In Eq. (33), the states are normalized as 

{B.,{v,\^)\B,{v,\^)) = 2«0(27r)3<53(k-k'), 

(i?;(z;,k,s)|B*(z;,k',s')) = 2v\2Trf5s,,5\\^-\^), 

(r,,(i;,k,s)|T,,(7;,k',s')) = v\2t,)H,s'5''{\^-\^), 

{S,,{v,\i,s)\S^j{v,'k',s')) = v''{2^f5ss'5''{\^-\^)- (35) 



B. Next-to-leading order analytic terms 



In this subsection, we investigate the NLO counterterms in the axial currents. Their matrix elements between single 
heavy hadron states are written as 



(Nf) 



in Eq. (32). These NLO counterterms play a significant role in the chiral expansion, since they have to be included 
to renormalise the one-loop contributions from the LO axial currents to matrix elements. 



First we notice that the chiral Lagrangian in Eq. \2Q\ does not contain any space-time derivative, therefore it does 
not lead to new terms in the axial currents upon applying the Noether theorem. To obtain the NLO axial currents, 
we introduce additional operators in the chiral Lagrangian, 



C 



(NLO, axial) 
HHxPT 



Ki^^ trD 
f K^^^ tro 



H 



kj 



Hf^ri, Hf\A,r\ str(xc) + nf^ tro [Rl'^n, hI'\A,U xf 
^if^ ep..p(5^z;''^"5^)f str (x?) + «f' e^-pl'^''^''^'') ^ str (A^x?) 



tro 



H?h'l,Hfx^{A^^ 



{S) 



AT) 



[{fA^S^)^ + {S^A^T)^] str(xe) 



(36) 
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where xs, is defined in Eq. (27). The mesonic sector of the above Lagrangian was already introduced in Refs. p51l27j . 
Upon applying the Noether theorem to Eq. (36), one obtains the currents which lead to the NLO analytic terms 

J^^h-^Ih- in Eq. (S, 



J.. 



(NLO, analytic) 



K""^ trD 



kl 



+ K^^-* trD 



kl 



' Kg tro 



kl 



str(x5) + nf^ trD 



<'7.7a Hf^ X,) 



^3 ^fii^o'P 



■Tl+lS'')^StTix^)+K['h^,,p{S''v'^S'')^Stv{Tl+l X?) 



AT) 



13, i Ati^l^)i 

(+)c ^ , /'e J+) 



{TXiAtls,)^ + {S,XirlllT)'^ 



{fr^lS,), + {S,rll^lT)^ str(xc) 



(37) 



where t^^^ is defined in Eq. (30). Although it is not explicitly shown in the above equation, these NLO currents 
depend on Nf. 



Comparing the currents Jij^l^'~'' ^^^^^'^^'^'^ to their leading-order counterparts, ^fj^"* in Eq. (29), one observes that they 
share the similar feature in the combination of the heavy-hadron fields with the flavour matrices rj^^^ . The complication 
in j^j^LO, analytic) j.gg^j^g Completely from the insertion of x^, which contains one power of the quark-mass matrix. 
This shows that one can write the NLO matrix elements as 

.(NLO, analytic) _ / W . I I ff. \ . „ ^ A/-(^/ ) 



NLO 



and 



O(TOq) - O(Mgoidstonc): 



where niq is the light-quark mass. 



IV. ONE-LOOP CONTRIBUTIONS IN SU(2) HHxPT 



We now turn to the discussion of the one-loop results for the axial-current matrix elements. In this section, we first 
present a simple case, namely SU(2) xPT in the infinite-volume limit, and use it to illustrate the main features of 
these one-loop contributions. Details of the SU(4|2) and SU(6|3) PQxPT results are addressed in the next section. 



We start by reducing the leading-order matrix elements in Eq. ( 33 ) to a simpler form. Notice that all these matrix 
elements are proportional to the axial couplings, (71,2.3 • Therefore, from the generic form of the chiral expansion for 
the axial-current matrix elements given in Eq. ( 32 ) , we can define the "effective" axial couplings 



(51) eff 

(52) off 

(53) cff 



51 X 



= 52 X 



.93 X 



1 / ^(2) ^(2) ^(2) p(2) \ ^(2) 



1 



f2 y ud 



(2) 



(2) 

Sdd 



p(2) 



n(2) 



(38) 



with the wavefunction renormalisation (W), tadpole (T) and sunset (Q) diagram contributions from Fig. [T](a), (b) 
and (c). 
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The result for the tadpole diagram is particularly simple. In the infinite-volume limit, it is 



following the definition of the function /(m) in Eq. (B4) in Appendix [B} Here is the pion mass, and fi is the 
renormalisation scale. The dependence on /i is cance led by the NLO counterterm contributions A/"^^^ in the above 
expression for the effective axial couplings. 

In this SU(2) full QCD case, the infinite-volume limit of the wavefunction renormalisation and sunset diagrams can 
be written in two functions 

Him, A) = 

Ai — A2 



with the function F defined in Eq. ( B4 ) . The scale A in these functions results from the mass difference between the 
external and the internal heavy hadrons. In the heavy quark and the isospin limits, we have 

Mb2 - Mb^ = Ms,, - Ms,^ = 0. 

Therefore the only relevant heavy-hadron mass difference in these limits is 

A(^) = Ms,„ - Mt,„ = Ms,, - Mt,^ ~ 200 MeV, 



and the effective couplings in Eq. ( 38 1 are 

(5i)eff = 31 1 



(52) cff 

(53) eff 



92 



93 



^I{M^) + j^H{M^,0) + analytic terms 



1 



2 

H{M^, 0) + ^ [h{M^, -A(^)) - 2K{M^, -A'-^\0)^ + analytic terms 
-2H{M^,A^^^)+H{M^,Q) 



+ ^ -A^-^)) + 9H{M^, A^^^) ~ 2K{M„, A^^\Q)^ + analytic terms 



(40) 



with the analytic terms resulting from M^"^^ in Eq. (38 1. Here we stress that the tadpole diagram is the dominant 
one-loop contribution to the chiral expansion of {gi)es- This is because the typical value of the coupling, gl ~ 0.25, is 
small, leading to the suppression of other diagrams in the above equation'^. A numerical comparison of the individual 



contributions from different types of Feynman diagrams will be given in Sec. IV B 



Before proceeding with further discussion of the formulae in Eq. (40), we notice that the function H{m,A) can be 
related to /(m) when A = 0, 



H{m,0) = -/(m) 



9 / ? 
m I m 

lovr^ 



This leads to the simplification of the chiral expansion of (gi)cff, 



(5i)off = 51 



;Mi log 



M'i 



(47r/) 



log 



Mi 



c{y)Ml 



(41) 



(42) 



^ Since 92,3 ~ this suppression is not present in the chiral expansion of (a2,3)efl- 
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The renormalisation-scale dependence from the loop diagrams is cancelled by the coefficient, c(/i), of the analytic 
term which also encodes the contributions from the NLO Lagrangian. 

In the following two subsections, we first address an issue related to the chiral limit of the formulae presented above, 
and then present an estimation for the numerical size of the one-loop corrections. 



A. Wavefunction renormalisation and sunset diagrams in the chiral limit 



As pointed out in Eqs. ( 39 ) and ( 40 1 , the infinite- volume one-loop contributions from the wavefunction renormalisation 



and sunset diagrams can 



the function defined in Eq. (|B4|), 
-1 



F{m,A) 



3e written in terms of the functions H and K, which are obtained by taking derivatives of 



167r2 



Alog ( ^ 



/10A2 



4m2\ ^ 2(A2-m2) /A\ 



where 



R{x) = \/ — I log (^x — \/x^ — 1 + iej — log (^x + \/ x^ — I + i 



This function is obtained by regularising the loop integrals with the subtraction scheme defined in Eq. (Bl) in 
Appendix |b] Implementing this scheme is a common practice in xPT calculations [28] . It leads to the result that 
F(m, A) does not vanish in the limit m unless A = 0. Such behaviour does not cause any conceptual problem in 
the effective theory, since the axial couplings, 51,2,3, can undergo finite renormalisation depending on the subtraction 
scheme used to regulate one-loop integrals. Various subtraction schemes always lead to the same physical quantities, 
such as the hadronic masses and axial transition amplitudes, which are scheme-independent. On the other hand, it 
would be desirable and natural to choose a scheme in which the one- loop contributions decouple in the chiral limit. 
As pointed out in Refs. [111110], it is possible to find a scheme such that the real part of F vanishes in the chiral limit. 
It is implemented by simply rewriting F as 



i^(""'^)(m,A) 



-1 



167r2 



m^Alog ( 



2A3 , / m2 

v°n4A^ 



4m2A 2(A^ 



-to2n /A\ 

'-mR - 

3 V TO / 



(43) 



and appropriately modifying the counterterms to absorb the difference (a finite polynomial in A). It is straightforward 
to demonstrate that when A + to > 0, in which case the external heavy hadrons are stable particles, this function is 
real and 



lim i^(™'^)(TO,A) ==0. 

m— ^0 



(44) 



In the case A + m < which corresponds to the situation that the external heavy hadron becomes unstable, the 
functions F and F'^''"'^^ are complex. Although the real part of F^^^^^ vanishes in the chiral limit, the imaginary part 
remains non-zero. This occurs when 



AU < \Mt- . - Mg^', I = A^-^) - 200 MeV. 



(45) 



Below this threshold, one cannot define matrix elements containing external S^, hadrons. In principle, more compli- 



cated matrix elements can be used to determine the couplings 172 and 173 for the pion masses in the regime of Eq. (451, 
but this is beyond the scope of this work. See Refs. [31] [32] for related discussions. Here we stress that one can 
perform lattice calculations in the regime where the pion mass is larger than A'^' but small compared to the chiral 
symmetry breaking scale, such that the external hadrons are all stable and the chiral expansion is still valid. These 
calculations enable the extraction the axial couplings, (71.2.3, which can then be used to perform chiral extrapolations 
and make predictions for other quantities. 



B. Evaluation of individual contributions 



In this subsection, we use the simple infinite volume, SU(2) case to explore the typical size of the one-loop contributions. 
This can be best summarised by the plots in Fig. [2] In these plots, the pion mass dependence of the loop contributions 
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FIG. 2; Comparison of the individual infinite-volume one-loop contributions to the pion mass dependence of the (real part of 
the) various effective couplings, {gi,2,3)^g evaluated using the values of the LECs given in the text. The kinks in the wavefunction 
renormalisation and sunset contributions to the baryonic couplings arise from the 5* — ^ Ttt threshold at = A'-^^ Below 



this threshold, the curves lose their physical interpretation. The subtraction scheme is that presented in Eq. ( 43 1 
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to three effective axial couplings [their real part in the case of {g2.3)^ff] is shown for exemplary values of the various 



IV A 



The 



low energy constants. These results are obtained using the subtraction scheme defined in Eq. (43) in Sec.^ 

leading order contribution is also shown. We take gi — 0.5, a value consistent with recent determinations |12H17j and 
then use the quark model expectations for the other couplings, g2 — 2gi and 53 — \/2gi (in our normalisation) [3]^ 
which are far less constrained. We work in the heavy-quark limit so that A^*^^ — 0, and we have set the S ~ T mass 
differences to A'^^ — 200 MeV, consistent with experiment [35]. The renormalization scale used here is fj, = 47r/. 

It is clear from these figures that the tadpole contributions provide an important part of the chiral non-analytic 
behaviour of the axial couplings. Furthermore, in the range of pion masses considered here, ^ 400 MeV, the NLO 
contributions from loops are numerically small corrections to the leading-order results.^ This indicates that in this 
range the SU (2) chiral expansion of the axial-current matrix elements is well-behaved. Variations of the low-energy 
constants, gi, 52, 53, A*^^\ and the renormalisation scale, within reasonable ranges do not substantially alter the 
behaviour shown in Fig. [2] 



V. ONE-LOOP CONTRIBUTIONS IN SU(4|2) AND SU(6|3) HHxPT 



In this section, we study the one-loop contributions in Eq. (32) in SU(4|2) and SU(6|3) partially quenched HHxPT in 
finite volume. These results are complicated because we keep the SU(3) light-flavour breaking effects from Eq. ( 26 ) in 
our calculation. Here we investigate the structure of the one-loop computation via analysing the quark flavour flow 
picture [31] . The details of the results are given in Appendices [C and [Dj 



A. The tadpole diagrams 



First we present the contributions from the tadpole diagrams. These take the simple form, 
-2I(M„,„0, 



-r(2) _ 

'ud — 



I I D O D 



where the functions I and I3 are defined in (B3) and (B13), respectively. The tadpole diagram results are completely 
determined by the structure of the axial currents. 



B. The self-energy diagrams 



In this subsection, we present the heavy hadron wavefunction renormalisation, resulting from the self-energy diagrams. 
These are more complicated than the tadpole-diagram results in our calculation, since we keep track of the flavour 



SU(3) breaking effects from both the Goldstone masses [Eq. ([3])] and the heavy-meson and baryon spectrum [Eq. (26)]. 
It is helpful to analyse the quark flavour flow diagrams [3j to understand the structure of the results. To investigate 
this structure, we first assign a "direction" to each flavour flow line: 

• The flow following the direction of a line means a quark with that flavour, while the flow against the direction 
means its anti-quark. 



* These values are also consistent with preliminary lattice QCD results 1331 . 

^ The imaginary parts of (g2,3)eff that arise for M,r < A^^^ are also small, | Im((ji2,3)effl < 0.05. 
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(a) 



(b) 



FIG. 3: Quark flavour flow structure for the meson seff-energy diagrams. The thick fine represents the anti-6 quark, and the 
thin fines are the vafence fight quarks, while the dashed fine is the sea fight quark. Diagram (a) contributes to the w and w' 
terms in Eq. (47 1 when the internaf heavy-fight meson is B*, and Bii respectively, while the Gofdstone meson is composed of 
a j vafence quark and an i' sea anti-quark. Diagram (b) is the "hairpin" structure and resufts in terms containing w (internaf 



B*) and w' (internaf Bj) in Eq. (47 1 



For the analysis of the heavy meson wavefunction renormahsation, we follow the nomenclature for the coefficients in 
front of the sum (integral) in a loop diagram: 

• The "tilded" coefficients accompany the "hairpin" contributions from the light flavour-singlet mesons. 

• The "primed" coefficients multiply the sums in which a B meson appears in the loop, while the "unprimed" 
coefficients are for the cases involving an internal B* meson. 

The quark flow picture for the heavy meson wavefunction renormahsation diagrams is presented in Fig. [3j Since there 
is only one valence light quark involved, and the internal valence-quark loops are cancelled by the ghost-quark loops, 
the only possible non-hairpin structure is from the sea-quark contributions. This is depicted in Fig. |3](a), where the 
Goldstone meson is composed of a j valence quark, and an i' sea anti-quark. The hairpin contribution is presented in 
Fig. 1(b). 

Following the above nomenclature and the quark flavour flow picture in Fig. [3j the results for the heavy-meson 
wavefunction renormahsation can be written as. 



(Nf) 



9l 



a -' 
a '- 

+ wi'^Jl niM,,a, 6^P) + *K -Hn, (M,,„, 0) , 



(47) 



where the summations are over the flavours u and u' in the SU(4|2) theory, and are over the flavours u, s, u' and s' 
in the SU(6|3) theory. The functions H and "H^V/ Siie results of the sums (integrals) involved in the loops, and are 
defined in Eqs. (B9), (Bll), \b12\ and (B13) in Appendixli] The mass parameters Mj^a, A^^^) and S^j^ are defined 



in Eq. (Al) in Appendix |A[ T'he coefficients w, w' , w and ur are presented in Table|l]in Appendix [C| 

Next, we discuss the structure of the baryon self-energy diagrams. We start by modifying the above rule for assigning 
the "primed" coefficients, 

• The "primed" coefficients multiply the sums in which the T baryon appears in the loop, while the "unprimed" 
coefficients are for the cases involving the internal S baryon. 

These diagrams are further complicated by the presence of two light valence quarks. To keep track of the flavour flow 
of these two quarks, we introduce an additional rule to our notation. 



• For a baryon (Tij or 5"^-,) of light flavour indices i and j, we assign the coefficient u to the diagram if the quark 
carrying flavour i appears in the Goldstone meson. For all the other cases, including the appearance of the 
anti-i in the Goldstone meson, they are accompanied by the coefficient w. 
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(a) 



(b) 



FIG. 4: Quark flavour flow structure for the baryon self-energy diagrams without the "hairpin" structure. The thick line 
represents the b quark, and the thin lines are the valence light quarks, while the dashed line is the sea light quark. Diagram 
(a) is the "crossing" type which does not involve sea quark contributions. The explicitly-shown diagrams give rise to the terms 
multiplied by w and w' in Eq. ( 48 1 when the internal baryons are Sai and Tai respectively [a = z in diagram (a) and a = z' in 



diagram (b)]. Interchanging the flavour indices i and j leads to the corresponding u and u' terms in the same equation. 






(a) 



(b) 



FIG. 5: Quark flavour flow structure for the baryon self-energy diagrams involving the "hairpin" structure. The thick line 
represents the b quark, and the thin lines are the valence light quarks. Diagram (a) contributes to the w (internal Sij baryon) 
and w' (internal Tij baryon) terms in Eq. (481, while Diagrams (b) and (c) result in the u (internal Sij baryon) and u (internal 
Tij baryon) terms in the same equation. 



The flavour flow structure for the baryon self-energy diagrams can be summarised in Figs. [4] and [5j For the diagrams 
exphcitly shown in Fig. |4] (a) and (b), the Goldstone mesons are composed of (j, anti-i) and (j, anti-i') respectively. 
Therefore, they are accompanied by the w-type coefhcients (w for the internal S baryon and w' for the internal T 
baryon). Terms with the u-type coefficients are obtained by exchanging the flavours i and j, as also indicated in this 
figure. Notice that the "non-hairpin" valence- valence Goldstone contributions appear via the "crossing" configuration 
in Fig|4](a). The "hairpin" structure of the baryon self-energy diagrams is presented in Fig.[5j From the above rules, 
it is clear that the diagram in Fig. [5] (a) leads to a w-type term, while those in Fig. [5] (b) and (c) are multiplied by 
w-type coefficients. 

Following the above discussion, we obtain the results for the baryon wavefunction renormalisation, 



{Nf 
T.J 



9l 



(Ni) 



9l 



E 



nM,,a, CI) + 4l'l ^N, [Aha. 0) 

+ ^(M...,-A(^^ +^^:^^ H^.,(M,„-A(-)) 



(48) 
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FIG. 6: Quark flavour flow structure for the meson sunset diagrams with external Bj and states. The thick line represents 
the anti-6 quark, and the thin hue is the valence light quark. Th e cross is the current Jjk- This diagram results in the y 
(internal Bj and B^,) and y' (internal B* and Bk) terms in Eq. (49 1. This is the only possible quark flavour flow conflguration 
for the meson sunset diagrams. 



where the summations are over the flavours u and u' in the SU(4|2) theory, and are over the flavours u, s, u' and 
s' in the SU(6|3) theory. The relevant coefficients, w, u . . . ar e presented in Tables [TT| and III in Appendix [C| The 
mass parameters Ma^b, ' - -'R"! , „ , . „ 

literature 



A^^^ and (JL L are defined in Eq. (Al) in Appendix A These results agree with those in the 



C. The sunset diagrams 



In this subsection, we discuss the structure of the sunset diagrams. The Lorentz indices carried by the hadronic states 
are completely absorbed into the tree- level contribution in Eq. ( 32 1 , therefore they are omitted in the notation below. 



In order to organise the results, we follow the same convention in assigning the "flow direction" to a quark line and 
the "tilded" coefficients to the terms involving the "hairpin" structure, as that in the self-energy diagrams. 

First we study the sunset diagram for the axial-current matrix element between the Bj and B'^ mesons. Because of 
the flavour-changing structure of the currents that we consider in this work, it is straightforward to demonstrate that 
in this case the Goldstone meson must involve the "hairpin" contribution. This is depicted in Fig. [6] Furthermore, 
the internal heavy meson with the light flavour j must be a B* since there is no _B—_B— Goldstone coupling in the 
Lagrangian or the current. On the other hand, the internal heavy meson involving the light flavour k can be either 
Bk or B^. These two cases are distinguished by the "primed" and the "unprimed" coefficients in the results. We then 
obtain the sunset-diagram contribution to this matrix element as 



9l 



-(Nf) 



IC^, (M,,,, A(^^), A(^^) + <5i^)) + y'l^j^. ic^, (M,,., A(^) , ^) 



(49) 



where /C and K-Nf are the sums (integrals) involved in the loops, and are defined in Eqs. (B9), (Bll ), (B12) and (B13) 
in Appendix [b| The mass parameters Mj^ki A*^*^) and 5^!'!^^ are defined in Eq. (Al) in Appendix |a| The coefficients 
y and y' are presented in Table VII in Appendix [D| 



Next, we investigate the sunset diagrams for the following axial-current transitions involving baryons. 



qfJ. 



St, Sr,, (50) 

where the spectator quark carries the flavour index i. The quark flavour flow configurations are shown in Figs. [7j 
[8] and |9] Again, we use the "tilded" coefficients to denote terms in which the "hairpin" structure appears. Because 
of parity, there are no axial couplings amongst even number of Goldstone mesons, therefore the flavour indices j 
and k must appear in the internal baryons. These internal baryons can be T or S* type. Denoting the other flavour 
index in the loop by a, we adopt the following convention to distinguish various possibilities for the internal S and T 
contributions: 



The SU(3) breaking effects arising from Eq. ( |26[ l are not included in the results in Ref. [35] . We have also checked these wavefunction 
renormalisation diagrams against the fuU-QCD, SU(3)-limit results at A'^' = in Ref. [5]) and found agreement. 
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(a) 



(b) 



FIG. 7; Quark flavour flow structure for the baryon sunset diagrams without the "hairpin" structure. The thick hue represents 
the b quark, and the thin hues are the valence hght quarks, while the dashed line is the sea light quark. The cross is the current 
Jjk- These diagrams lead to the x (internal Saj and Sak baryons with a = i' ,j, k in diagram (a), (b), (c) respectively), x' 



(internal Taj and Sak baryon ) and x" (internal Saj and Tak baryons) terms in Eq. (51 1 






(a) 



(b) 



(c) 



FIG. 8: Quark flavour flow structure for the baryon sunset diagrams with the "hairpin" structure involving the spectator quark 
(flavour i). The thick line represents the h quark, and the thin lines are the valence light quarks. The cross is the current Jjk- 
These diagrams lead to the x (internal Sij and Sit baryons), x' (internal Tij and Sik baryon ) and x" (internal Sij and Tik 
baryons) terms in Eq. 1 51 1. 



• If the internal baryons are Saj and Sak, then the coefficient for the diagram is "unprimed". 

• If the internal baryons are Taj (left in the loop) and Sak (right in the loop), then the coefficient for the diagram 
is "primed" . Such terms are absent in the T S transition amplitudes. 

• If the internal baryons are Saj (left in the loop) and Tak (right in the loop), then the coefficient for the diagram 
is "double-primed". 

To keep track of the flow of the spectator quark i in these processes, we follow the rules: 

• If the spectator quark flavour is present in the Goldstone meson, then the diagram corresponds to a term with 
x— type coefficient (x, a;', x", i, x' or x"). 

• If the spectator quark flavour is absent in the Goldstone meson, then the diagram corresponds to a term with 
?/— type coefficient. 



In Fig. [7j we show the quark flavour flow diagrams containing no "hairpin" structure. In such flow configurations, 
the spectator quark flavour always appears in the Goldstone meson. Therefore they will only be accompanied by 
the a;— type coefficients. Notice that the valence-valence Goldstone mesons also appear in these diagrams via the 
"crossing" configurations in Figs. [T] (b) and (c) . The "hairpin" contributions to the quark flavour flow configurations 
for the processes in Eq. (50) are shown in Figs, [s] and [o] These two figures are distinguished by the presence/absence 
of the spectator quark flavour in the Goldstone propagator. Therefore they correspond to terms with x— and y— type 
coefRcients respectively. 



Following the above rules in analysing the quark flavour flow structure, the results for the baryon sunset diagrams 
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FIG. 9: Quark flavour flow structure for the baryon sunset diagrams with the "hairpin" structure in which the spectator quark 
(flavour i) is absent. The thick hne represents the b quark, and the thin lines are the valence light quarks. The cross is the 
current Jjk- These diagrams lead to the y (internal Sij and Stk baryons), y (internal Tij and 5*;*: baryon ) and y" (internal 
Sij and Tik baryons) terms in Eq. (511. 



can be written as 





= 9l< 
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^ a 












+9l< 


^ a 


^"(Nf) 






: ,-/'(JV/) r 


q(Ni) 


= 9l< 
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^ a 


Sij, Sik -a 




+9l< 




Sij -Sik ,0- 












E 

a 


S ij : Si k 1 



(51) 

where the summations are over the flavours u and u' in the SU(4|2) theory, and are over the flavours u, s, u' and s' in 
the SU(6|3) theory. The relevant coefficients are presented in Tables [Tvl [v] [yi} [Vnl |VIII[ [DC} and|x]in Appendix |d] 



VI. Hi Hi tt{K) transition AMPLITUDES 



The axial-current matrix elements, presented in the previous sections, are closely related to those in the strong-decay 
amplitudes, such as 

B* ^ B„7r, B* ^ BuK, 



.(*) 



AfcTT, 



(52) 



Note that, with the exception of SJ, ^ A^tt, for bottom hadrons the above decays are kinematically forbidden in 
nature. In HHxPT, the LO and NLO analytic terms for these decay amplitudes have the same structure as the matrix 
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(a) 



(b) 



FIG. 10: Diagrams contributing to the decay amplitudes in Eq. 1 52 1. The self-energy diagrams leading to wavefunction 
renormalisation of the external particles are not shown in this figure. The dashed lines are the Goldstone mesons. The single 
solid lines denote generically the external heavy hadrons, while the double solid lines are the internal heavy hadrons. They can 
be B, B* mesons or Ti-j, Sa baryons. The vertices are all from the axial-coupling terms proportional to gi,2,3 in the strong 



chiral Lagrangian in Eq. (241. Diagrams (a) and (b) are the "tadpole" and "sunset" types respectively. 



elements in Eq. (32 1. That is, the LO contributions are all proportional to the axial couplings 51,2.3, while the NLO 



results are polynomials in the Goldstone masses. Therefore we only address the one-loop diagrams for these decays. 



To compute the one-loop amplitudes for the processes in Eq. ( |52[ ), one has to calculate the wavefunction renormali- 
sation of the Goldstone bosons and the heavy hadrons, as well as the tadpole and sunset diagrams in Fig. [TOj The 
Goldstone boson wavefunction renormalisation can be found in standard references such as [Bj and i22j , and the heavy 

^ (|48l 



hadron wavefunction renormalisation is presented in Eqs. (47) and (|48|. The amplitudes from the sunset diagram in 
Fig. 10 (b) are identical to those from the corresponding diagram in Fig.[l](c). Therefore they are equal to the results 



and (51 1. The tadpole diagram in Fig. 10 (a) differs from that of the axial-current matrix 



presented in Eqs. 



elements in Fig.[T](b) by a factor of one-third. That is, one can take the results in Eq. (46), and multiply them by 1/3 



to obtain the corresponding tadpole-diagram contributions to the decay amplitudes in Eq. (52). It turns out that the 



contribution from the tadpole diagram is exactly cancelled by the contribution from the wavefunction renormalization 
of the external Goldstone boson [35' . As is shown in Fig. [2] the tadpole diagrams provide significant contributions to 
the axial current matrix elements and will lead to significant differences between the quark-mass dependence of axial 
current matrix elements and that of strong decay amplitudes. 

These decay amplitudes have also been computed in Ref. to one-loop order in SU(3) HHxPT in the infinite- 
volume limit with A^*^) = A^-^^ = 0, and without the SU(3) breaking effects from the Lagrangian in Eq. (26). Our 
results agree with those presented in Ref. (36J in the same limits. 



VII. CONCLUSION 



With the expectation of precise data from the LHCb collaboration and from the potential SuperB experiment, accurate 
QCD calculations of quantities involving B mesons and single-6 baryons will be important in further constraining 
flavour physics and in looking for physics beyond the SM. This is a challenging but necessary task. In this paper, we 
have presented calculations for axial-current matrix elements involving single heavy hadron external states in HHxPT 
at the NLO. We have performed these computations in partially quenched xPT for both Nf — 2 and Nf =2 + 1, 
including finite volume effects. Our results are essential for extracting the axial couplings in HHxPT from experimental 
data or lattice QCD. These axial couplings are central quantities in b physics, as they control the light quark mass 
dependence of &-hadron observables and determine the strong decay widths of heavy hadrons. 

We have discussed the SU{2) case in detail, numerically analysing the behaviour of the various loop contributions 
for natural values of the low-energy constants. Based on our study, we conclude that the SU{2) chiral expansion of 
the axial current matrix elements is well-behaved for A/^r ^ 400 MeV. This implies that lattice calculations that are 
performed in this regime can be used to determine the axial couplings reliably. 
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Appendix A: Mass parameters 

In this appendix, we define various quantities appearing in our results. First we present the hadron masses. 

Ml, = Bo {ma +mt), ~S'ys = „ - Af^ „, , S^ss = - Ml,, , Mj- = ^ (m^ „ + 2MI, - 25ls " ^Kss) , 
A(^'^) = Afs* - Mb, , Sl'^p = Mb^ - Mb, = A/b* - Af^. , 

A(^) = Afs„. - Mt„, , <5if,L = Mt^. - Mt.. - Ms^. " ^5.. , (AI) 

where Bq is defined in Eq. m. As explained in the main text, A*^^^^ vanishes in the heavy quark limit, while A^^^ 
remains non-zero and is of 0(Aqcd)- In this paper, we work in the isospin limit, and denote the pion mass as Af„_„. 

It is useful to define the following quantities which appear in the "hairpin" contributions to the flavour-singlet meson 
propagators in the SU(6|3) theory. 

\ 2 

r2 



2 {S^vs - + Mi) 3 _ 3 (8~5t.ss + {2~Sh " + A/,,, 



^ 

_ 24, ^ ^ 6^,, {2~5h^Ml^ + Mt) 

~ M2 - Ml ' ~ op + 4^2 - M2 + M2 ' 

n(«) — ^ L n^"' — ^ ^ 



(Af2 „ - Ml) [Ml^ -Ml)' "-^ (Af 2 „ - Ml) {Ml - M^ j 



{mI.u - Ml - 2py,) (mi - Ml - 25lss. 

{Ml^~Ml){Ml-Ml) ■ ^ ^ 



In the full QCD limit, where m„' = m„ and 



^QCD ^ 3 -QCD ^ o (^QCD ^ n pQCD ^ n 

i^C-OQCD^O, i?(^)QCD^o, i^WQCD^^^ (A3) 
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Appendix B: Integrals and sums 



In this appendix, we present results of loop integrals and sums using dimensional regularisation, with the ultra-violet 
divergences removed by subtracting the term, 



A 



2 



-7£;+log(47r) + l, 



(Bl) 



where d is the number of space-time dimensions. This is a commonly- used scheme in xPT calculations [28^. It is 
different from the MS scheme by the constant "1" on the right-hand side of the above equation. It can also be changed 
into the scheme discussed in Sec. IV A straightforwardly. Finite volume effects in the limit mL ^ 1 (rn is a generic 



Goldstone mass and L is the spatial lattice volume) are computed by replacing the momentum integrals by sums in 
the spatial directions. The one-loop contributions appearing in this work can all be obtained by investigating the 
following sums /integrals 



I{m) 



,4-d 



d'^k 



TO 



J-(to,A) EE {gP'^-vPy") 



[d^)4^j27T 



ie IGtt^ 



A, 



(B2) 



{d~l) ^ {2ttY (fc2 - m'^ +ie){v ■ k - A + ie) 
where n is the renormalisation scale, and the symbol 



167r2' 



2A2 



t 



d^k 



means performing the sums in three spatial directions using the Poisson summation formula, followed by dimensionally 
regularising the infinite-volume integrals. 

We can further separate the infinite-volume limit of X and T from the finite- volume contributions, 

I(to) = I{ni) + Ipvim), 



T{m) ^ F(m,A) + FFvM- 
The functions / and F are results from the ordinary one-loop integrals. 



I{m) 
F{m,A) 



w? ( rr? 



167r2 



2 2A' \ ^ , /TO- 

Ulogf-, 



/10A2 4to2 
V^9 3^ 



with 



, 2(A2-to2) /a 

A -I- -mR — 

3 V TO 



Rix) = \/ ~\ log [x — \/x2 — 1 + it^ — log + \/a;2 — 1 + ie^ 



(B3) 



(B4) 



(B5) 



The funct ion _F(to, A) does not vanish in the to — > li mit u nless A = 0. One can adopt the scheme discussed in 
Sec. IV A by simply rewriting F as defined in Eq. (43), and the real part of the function F^^^^'^ is zero in the 



chiral limit for arbitrary A. 
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For the case in which the external hadrons are stable particles, the finite- volume pieces can be shown to be [SJ [57] 

m 



hvim) 



A 9 ~T^i [umL) 
47r2 ^ uL ^ ' 



127r^ uL Jq 



X <^ 1 + 



|fc|sin(u|fc|L) 



3 15 
8umL ^ 128{umLY 



O 



1 

umL 



^2 ^~umL 
iL>l —TO \ - e 



A 



E 



247r ^ uL 



A. 



where u = (711,1*2,^3) with itj G Z, u = |m| and 
= e("')[l-Erf(z) ] 



umL 



( ' ' 








\umL ^ 


) 






39z 






"1^ 
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2z2 ) e(^')[l-Erf(z) ] 



(B6) 



(B7) 



-^ + ^)e(^^)[l-Erf(z)] 



with 



A\ I umL 

m , 



(B8) 



Higher order terms in the \/{umL) expansion in Eq. (B7) can be easily calculated. The integer Ui can be interpreted 
as the number of times that the pion wraps around the spatial volume in the i— the direction. 

The functions appearing in our one-loop results are 



dT{m, A) A A ^ -F(m, Ai) - J-(m, A2) 

n{m, A) = and /C(m, Ai, A2) = — 



Ai- A2 



and /C^'(m, Ai, A2) = 



_ a/C(m,Ai,A2) 



dm^ 



Notice that 



JC (m. A, A) EE ^Ihn^ JC{m, A, A') = H (m. A) . 



(B9) 



(BIO) 



To present the residual flavour-singlet "hairpin" contributions in a compact form, we define three functions 

iNf{m), nNf{m,A), and /Cjv^.(to, Ai, A2). (Bll) 

They take the explicit form 

12 (m) = I{m) + 2Sys Ijj'{m), 



H2{m,A) = n{m,A) + 2Sls'H^'{m,A), 
^2(to,Ai,A2) = /C(m,Ai,A2) -t- 252,5 /C^,(m,Ai,A2), 



(B12) 



Similar formulae for finite- volume effects are also obtained in Ref. |38| . 
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TABLE I: Coefficients for heavy-fight meson wavefunction renormafisation, in Eqs. (47 1, in the isospin limit 



in the SU(4|2) theory (iV/ = 2) where the mass m in the arguments is always equal to M„_„, and 
+ (1 - 5,.,) [d^I I(A/„,,) + D^^ I{M,.,) + I{Mx)] , 

n3iMa,b,A) - Sa.b{Aa,bniMaM,^) + il-Aa,b) H ( A/x , A) + ^.fe ( Af„,b , A) } 

+ (1 - 6,,,) [dI""} n{Ma,a,A) + dI'I n{M,^b, A) + Dif n{Mx,A)} , 

^3(Ma,6,Ai,A2) = Sa^i,{Aa,blC{Ma,b,Al,A2) + il-AaM) /C(Afx , Ai , A2) + Ca,6 /C^' ( Af,^6, Ai , A2) } 

+ (l-5a,6){^5 /C(A/4^„,Ai,A2)+i^i'i/C(A4,,,Ai,A2) + i?2' ^(Afx, Ai, A2)} , (B13) 
in the SU(6|3) theory {Nf = 3). 

Appendix C: Coefficients for wavefunction renormalisation 



In this Appendix, we present the coefficients in Eqs. (47) and (48 1 relevant to the matrix elements investigated in 



this work. These coefficients are summarised in Tables |I[ |II| and HI Because of isospin symmetry, it is not possible 
(or necessary) to distinguish between the w and u coefficients for some of the hadrons in the current study. For such 
cases, we simply present -u; + m in the tables. 
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nts for Tij baryon wavefunction renormalisation, 
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TABLE III; Coefficients for Sij baryon wavefunction renormalisation, in Eq. (48 1, in the isospin limit 



Appendix D: Coefficients for tlie sunset diagrams 



In this Appendix, we present the coefficients in Eqs. (491 and (51 1 relevant to the matrix elements investigated in 
this work. Because of the isospin symmetry, it is impossible to distinguish between some y coefficients and their x 
counterparts. For such cases, we put the symbol +x in Table VII and then present x + y in Tables VIII[ IX and |X] 
in the form that x is written as (number - y). 
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TABLE IV: Coefficients x\j in Eq. |5ll, in the isospin limit. 
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TABLE V: Coefficients 
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in Eq. 
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isospin 


limit. 
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TABLE VI: Coefficients 


'(^/) 


in Eq. 





, in the 


isospin 


limit. 
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TABLE VII: Coefficients y^jj/lj^, VhI^h^ ^^"^ fe/^H^i Eqs. (49l and JsiL Due to the isospin symmetry, some of these 

(Nf) 

coefficients cannot be distinguished from their H2 counterparts. For such cases, the vaiues are denoted +i in the table, 
and are presented together with the corresponding x^h/h^ ■ 
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TABLE VIII; Coefficients x^_^ in Eq. Jsik. Due to the isospin symmetry, some of them cannot be distinguished from their 

(N f) -.(Nf) (Nf) 

yj{{H2 counterparts. For such cases, we present Xg_^ + ijjj^ jj^ in the table. 
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TABLE IX: Coefficients ^^^^ ^"-l- ^^If- isospin symmetry, some of them cannot be distinguished from their 

Vhi H2 counterparts. For such cases, we present x^_^ + y^_^ in the table. 
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TABLE X: CoefBcients in Eq. |5l|. Due to the isospin symmetry, some of them cannot be distinguished from their 

-I'iNf) _^ -"{JVf) -"(Nf) . . . , 

H2 counterparts, tor such cases, we present + y^_^ m the table. 
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